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I.  INTRODUCTION 


Quantitative  reconstruction  of  climatic 
parameters  from  both  oceanic  and  continental 
records  shows  several  time  scales  of  the  suc¬ 
cessive  glacial-interglacial  episodes  that  have 
characterized  the  Earth  s  climate  about  the 
past  million  years.  For  example,  after  analyz¬ 
ing  data  from  the  2,083  m  ice  core  recovered 
by  the  Soviet  Antarctic  Expeditions  at  Vostok 
(East  Antarctic).  Jouzcl  ct  al.  (1982)  found  se¬ 
veral  peaks  (25.3,  45.7,  and  107.5  kyr)  in  the 
variation  spectra  of  the  Vostok  isotope  tem¬ 
perature,  and  related  these  peaks  to  the  astro¬ 
nomical  forcing,  i.e.,  the  obliquity  of  the 
Earth's  axis  (period  at  41  kyr)  and  the 
precession  of  the  equinox  (periods  at  23  and 
19  kyr).  Is  the  climatic  variation  purely  caused 
by  the  external  forcing?  If  so.  the  climate  pre¬ 
diction  would  be  relatively  simple  because  it 
becomes  the  forecast  of  the  Earth  s  orbit  tilt 
and  precession.  In  fact,  the  climate  system  is 
much  more  complicated.  The  climatic  variation 
is  caused  by  both  external  and  internal  sources 
( Fig.  1). 
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Fig.  1  A  simple  model  of  the  climate  system  il¬ 
lustrating  internal  and  external  sources  of  cli¬ 
matic  variation  (after  Imbrie,  1982). 
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The  internal  forcings  arc  processes  and 
mechanisms  within  the  Earth  s  system,  giving 
rise  to  self-excited  climate  fluctuations.  The 
temperature-albedo  feedback  mechanism  asso¬ 
ciated  with  long  and  short  wave  radiation 
fluxes  at  the  Earth's  surface  arc  given  consid¬ 
erable  attention  in  glacial  and  climate  models 
(e.g..  Kallen  et  al.,  1979;  Oerlemans  and  van 
der  Veen,  1984).  However,  not  until  recentlv 
(Chu,  1990).  the  role  of  hydrological  cycle  on 
the  climatic  variation  on  glacial  and 
interglacial  time-scales  is  given  less  attention 
although  it  is  realized  that  ice  freeze  melt  and 
evaporation/precipitation  are  contributors  to 
the  change  of  air  and  ocean  temperatures. 
Taking  hydrological  cycle  into  consideration, 
several  intrinsic  time  scales  of  the  atmosphere- 
ocean-cryosphere  climate  system  will  be  dis¬ 
cussed  in  this  paper. 

2.  ATMOSPHERE-OCEAN-CRYOSPHERE 
CLIMATE  SYSTEM 

The  total  surface,  ice  covered,  and  open 
ocean  areas  are  set  to  be  1,  N„  and  1  -  re¬ 
spectively  (Fig.2).  The  interactions  among  the 
atmosphere,  ocean,  and  ice  is  proposed  by  Chu 
(1990)  as  follows.  From  the  global  point  of 
view,  ice  advance  isolates  the  ocean  from  the 
atmosphere  and  reduces  the  heat  loss  from  the 
ocean  surface,  which  in  turn  decreases  the  air 
temperature  and  increases  the  ocean  temper¬ 
ature.  The  air  and  ocean  feed  back  into  the  ice 
through  two  different  kinds  of  mechanisms 
(Fig.3):  (1)  thermodynamical  feedback;  the  in- 
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crease  of  ocean  temperature  melts  the  ice  and 
causes  the  ice  retreat  (negative  feedback  be¬ 
tween  ice  and  ocean),  however,  the  decrease  of 
air  temperature  causes  freezing  and  makes  the 
ice  further  advance  (positive  feedback  between 
ice  and  air);  (2)  mechanical  feedback;  the 
modification  of  the  atmospheric  and  oceanic 
temperatures  varies  the  surface  evaporation 
rate,  which  in  turn  changes  the  sea-level  height 

and  the  ice  accumulation  rate.  Both  effects 
lead  to  a  further  change  of  the  ice  flow. 


The  atmosphere  is  heated  (or  cooled)  by 
radiation,  and  heat  fluxes  (including  latent  heat 
flux)  from  the  air-ocean  and  air-ice  interfaces, 


where  Tc  is  the  air  temperature,  is  the  air 
specific  heat,  Qat,  are,  respectively,  the  heat 
fluxes  across  air-ice  and  air-ocean  interfaces, 
Ra  is  the  net  radiation  absorbed  by  the  atmos¬ 
phere,  and  p„  and  H0  are  the  characteristic  val¬ 
ues  of  atmospheric  density  and  thickness. 

In  the  present  study,  it  is  assumed  that  the 
water  in  the  atmosphere  is  always  in  the  vapor 
form.  As  soon  as  condensation  happens,  the 
condensed  water  droplets  arc  assumed  to  be 
precipitation  right  away.  Therefore,  the  time 
rate  change  of  specific  humidity  q,  is  caused  by 
the  excessive  surface  evaporation  over  precipi¬ 
tation  and  ice  accumulation 
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Fig. 2  Physical  processes  in  the  coupled  atmos¬ 
phere,  ocean,  and  cryosphere  climate  system 
(after  Chu,  1990). 
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where  E,  Pr,  and  cs  are  rates  of  ocean  surface 
evaporation,  precipitation,  and  ice  accumu¬ 
lation.  pw  and  p ,  are  the  characteristic  values 
of  ocean  and  ice  densities. 
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2.2  Ocean 
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Similar  to  the  atmosphere,  the  ocean 
thermodynamics  can  be  depicted  by 
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Fig. 3  Feedback  mechanisms  among  air,  ice 
and  ocean. 
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where  T.  is  the  ocean  temperature,  c „  is  the 
ocean  specific  heat,  is  the  heat  fluxes  across 
ice-ocean  interfaces.  R.  is  the  net  radiation  ab¬ 
sorbed  by  the  ocean,  and  H.  is  the  ocean 
thickness. 

From  a  thermodynamical  point  of  view, 
the  ocean  thickness  is  reduced  by  the  excessive 
surface  evaporation  rate  £  over  precipitation 
rate  Pr,  and  by  the  ice  edge  freezing  rate  £. 

-|j-[(l-  SdpJU  -  ~  (»  -  Ai )PJE~  Pr) 

-  ( 1  -  N,)pA(F)F -  NtPA(  -F)F  (3) 

where  A  is  the  Heaviside  function,  i.e., 
A(*)-l,  as  xSO,  and  A(x)  =  0,  as  x<0. 
£>0  .is  taken  as  ice  freezing;  and  £<0.  is 
taken  as  ice  melting. 

2.3  Ice- shelf 

The  ice  is  assumed  to  be  isotropic  with 
thickness  II,  (Fig. 2),  and  the  ice-spreading  rate 
is  computed  by  (Wcertman.  1957) 

where  //,  is  the  height  between  the  ice  top  and 
the  sea-level  (Fig.2).  Values  of  n  vary  from 
about  1.5  to  4.2  with  a  mean  of  about  3,  and 
for  randomly  oriented  polycrystalline  ice  at 
-  |0*C  and  n  =  3.  a  value 
t4  =  3  x  \Q-*yr-'kPcr'  is  reasonable  for  ice 
stress.  The  angle  brackets  denote  values  aver¬ 
aged  over  ice  thickness. 

The  change  of  ice  coverage  .V.  is  due  to 
ice-spreading  and  ice  freezing  melting: 


A,  <£„>  + 


(I  -  S,)\(F\F 

H, 


NA(  -F)F 


(5) 


Following  Killworth  (1979).  the  heat  flux  is 
accomplished  through  an  exchange  coefficient 
K  and  is  of  the  form  K{  T  -  7}).  where  T,  re¬ 
presents  the  freezing  point  of  the  sea  water. 

Ihc  freezing  rate  is  therefore  computed  by 

T$+pj*KJTm-  7})] 

(6) 

where  K. .  Kg  are  the  heat  exchange  coefficients 
for  water-ice  and  air-ice.  respectively.  L,  is  the 
latent  heat  of  ice.  If  T,.  T,  arc  greater  than  Ti. 

F  is  taken  as  a  negative  value,  which  indicates 
ice  melting. 

Based  on  the  continuity,  Shumskiy  (1965) 
proposed  a  method  to  determine  the  change  of 
ice  thickness  by  comparing  the  accumulation 
rate  with  horizontal  strain-rates.  Neglecting 
nonlinear  inertial  terms  the  time  rate  change 
of  ice  thickness  becomes 

dH, 

~  cs  ~  Hi  <  tf,  >  (7) 

The  atmosphere-ocean-cryosphcrc  climate  sys¬ 
tem  is  depicted  by  ( 1 ),  (2).  (3).  (5),  and  (7). 

2.4  Interfacial  Heat  Fluxes 


The  heat  fluxes  across  the  air-ocean,  air- 
ice,  ice-ocean  interfaces  change  the  thermal 
features  of  each  component  of  the 
atmosphere-ocean-  cryosphere  climate  system. 
There  are  many  different  methods  to  compute 
these  fluxes,  among  them,  the  simplest  form 
was  proposed  by  Killworth  (1979): 

Qa  =  P^paFJ  T,  -  Tt) 

Qtra  =  PtfpaKJJ*  -  Fa) 

Qwi  =  PJCpwKU  Tw  -  Ti)  (8) 
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3.  ROLE  OF  THE  HYDROLOGICAL  CY¬ 
CLE 

Summation  of  (lb),  (3),  (5)  multiplied  by 
H„  and  (7)  multiplied  by  N,  leads  to 

N,)Pwllw  +  pMN,  +  J  =  0(9) 

which  shows  the  conservation  of  the  water 
mass  in  the  Earth  system. 

The  water  vapor  evaporated  from  the 
global  ocean  can  be  estimated  by 

E=CD\V.\j^UiTw)-q,-\  (10) 

where  CD  is  the  drag  coefficient,  V„  is  the  sur¬ 
face  wind,  q,  is  the  saturated  mixing  ratio,  and 
q,  is  the  mixing  ratio  near  the  ocean  surface. 
If  the  Earth  surface  is  assumed  to  be  covered 
by  ocean  and  ice  only,  the  water  vapor  evapo¬ 
rated  from  the  ocean  surface  is  transported  by 
the  atmospheric  circulation,  condensed  and 
precipitated.  The  precipitation  over  ice  leads  to 
ice  accumulation,  therefore,  if  the  total  water 
vapor  evaporated  from  the  ocean  surface  is 
balanced  by  the  precipitation  over  ocean  (Rr) 
and  the  ice  accumulation  rate  over  ice  (cr),  the 
relationship  among  surface  evaporation  rate, 
precipitation  rate  (over  ocean),  and  ice  accu¬ 
mulation  rate  should  be 

(\-Nl)pw(E-Pr)  =  N,pics  (11) 


discussed  in  this  paper,  the  external  forcings 
mainly  change  the  radiation  absorbed  by  at¬ 
mosphere  (R„)  and  ocean  (RJ. 

On  the  other  hand,  the  internal  forcings 
are  processes  and  mechanisms  within  the  Earth 
system,  giving  rise  to  self-excited  climatic  fluc¬ 
tuations.  Substitution  of  (8)  into  (1)  and  (2), 
leads  to  thermal  relaxation  lime  scales  for  the 
atmosphere 

+  (12) 

and  for  the  ocean 
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The  ice-spreading  time  scale  can  be  defined 
from  (5), 

TJ*=(<Cft>rl  (i4) 

Horizontal  ice  growth  (or  ice  edge 
freezing: melting)  time  scale  is  obtained  from 
(5): 


ni 

f  F[A(F)-A(-FK 

The  perturbed  ice  accumulation  c\  is  caused 
by  the  perturbed  evaporation,  precipitation, 
and  the  ice  advancement: 


4.  TIME  SCALES  OF  THE  CLIMATE  SYS¬ 
TEM 


(1  ~  A !,)PW 
N,Pi 


( E’-Pr')-lpw(E-Pr)  +  plcsl\ ' 


In  explaining  the  climatic  variability,  one 
may  first  ask:  is  climatic  variability  caused  by 
external  or  internal  forcings?  The  external 
forcings  are  agencies  outside  the  Earth  system, 
such  as  astronomically  caused  changes  in  the 
Earth's  orbital  parameters.  In  the 
atmosphere-ocean-cryosphere  climate  system 


(15) 

The  perturbations  of  evaporation  and  precipi¬ 
tation  rates  E'  and  Pr'  van-  on  a  much  shorter 
time  scale  than  the  tiine-scale  for  ice 
advance  retreat.  Another  time  scale  for  the  ice 
growth  in  vertical  direction  can  be  defined  by 

P(H, 

!>..(£ -Pr)  +  p,cJ  (l6) 
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Here  the  variables  without  means  the  basic 
states.  By  the  use  of  the  values  of  the  param¬ 
eters  listed  in  Table  l,  these  time  scales  arc  es¬ 
timated  as: 


the  low  frequency  modes  will  be  generated.  If 
the  coupled  system  perturbed  from  its  equilib¬ 
rium  states,  after  mathematical  manipulation, 
the  perturbation  equations  arc 


ta  =  1 1.6*/,  tw  i  66  yr. 

r,  ~  400it,  Tj  ~  3000>t 

Since  the  atmospheric  and  oceanic  thermal  re¬ 
laxation  times  scales  arc  so  much  shorter  than 
the  ice-spreading  time  scale,  the  atmospheric 
and  oceanic  temperature  perturbations  T.  and 
r.  almost  "instantaneously"  follow  the  ice 
advance  retreat  processes  in  the  aimosphcrc- 
ocean-cryosphere  climate  system.  From  ( 1 )  and 
(2).  we  may  obtain 

Tt-~  A„raiyif  =  A.r^V,  ( 1 7) 


where 

A.  =  (£U  4  Aw  = 

Equation  (17)  clearly  shows  the  cooling  of  the 
atmosphere  ( T,  <  0)  and  warming  of  the  ocean 
(7\  >0)  due  to  the  ice  advance.  The  per¬ 
turbed  ice  freezing  rate  becomes 

F  -  -  -  a'A.TjY,  (18) 


where 

=  Pacpa^iJll^Pl^l^l)’  =  Pw^pw^wlKPl^'J 


.  The  perturbed  ice  spreading  rate  is  computed 
by 


H: 


(19) 


-I  h'l 

-MlA'it (20) 

t-^+d+^r’]-^- 

-0.l(.1rr,-A7'tr')-V,  (21) 


where 


PJL-PM 

^ -X)  “2J3‘ 


and 


The  general  solutions  of  equations  (20) 
and  (21)  have  the  following  forms: 

<*7  (22) 


Substitution  of  (22)  into  (20)  and  (21)  leads  to 
a  second  order  algebraic  equation  for  o> 

u2  4-  io>[(  1  -  a)rj]  4  (1  +  fi2  4  A>,m3)t7’]  - 


6.  LOW  FREQUENCY  MODES  GENER¬ 
ATED  BY  THE  FEEDBACK  MECHANISM 


[(1  -er)(l  4m2)t71t/'  4  nxit3Nfa2  4  “  0 


All  the  lime  scales  presented  in  the  last 
section  arc  shorter  than  the  time  scales  of  the 
successive  glacial-interglacial  episodes  (10-100 
kyr)  shown  in  Fig.l.  If  we  consider  the 
atmosphcre-occan-cryos  phere  as  a  one  sy  stem, 


(23) 

The  frequency  here  is  a  complex  number,  i.c., 
cj  =  cj.  4  ioi„  to,  is  a  growth  rate,  and  u,  is  the 
frequency. 


437 


Period  (kyr) 


I 


7.  FREQUENCIES  FOR  NEITRAL  STA¬ 
BILITY  MODES 

The  modes  with  zero  growth  rate  (10  =  0) 
arc  called  neutrally  stable  modes.  The  solution 
of  (23 )  with  at,  =  0  is 

w,  =  -  [(1  +  ji2)3  +  A>,/i2^3]}l/t24) 

where 


The  time  scale  for  the  atmosphere-ocean- 
cryosphere  climate  system  depends  on  the  ratio 
of  the  two  shorter  time  scales:  r,  and  r,.  Fig.4 
shows  the  dependence  of  periods  (2^/cd.)  of  the 
neutrally  stable  modes  on  the  ratio  r.  When 
v  —  7.2.  the  period  is  43.3  kyr.  which  indicates 
that  the  feedback  mechanisms  among  atmos¬ 
phere.  ocean,  and  cryosphcrc  can  also  generate 
climate  change  on  10- 100  kyr  time  scales. 


Fig.4  Periods  versus  r,/r,  for  the  marginal  sta¬ 
ble  modes. 


DISCUSSION 

The  conventional  theories  refer  the  climate 
change  on  10-100  kyr  time  scales  to  the  con¬ 
centration  of  variance  near  the  Earth  s  orbit 
tilt  and  precession  frequencies  (Jouzcl  ct  al., 
1987).  However,  the  conceptual  atmosphcrc- 
ocean-cn  ospherc  model  described  in  this  paper 
shows  a  possible  positive  negative  feedback 
mechanism,  induced  by  the  hydrological  cycle. 
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